RUNHETC-2006-13 



Ising Spectroscopy I: Mesons at T < T c 

P. Fonseca and A. Zamolodchikov 1,2 

^^NHETC, Department of Physics and Astronomy 
Rutgers University 
Piscataway, NJ 08855-0849, USA 

2 L.D. Landau Institute for Theoretical Physics 
Chernogolovka, 142432, Russia 



Abstract 

This paper is our progress report on the project "Ising spectroscopy" , de- 
voted to a systematic study of the mass spectrum of particles in the 2D 
Ising Field Theory in a magnetic field. Here we address the low-temperature 
regime, and develop a quantitative approach based on the idea (originally 
due to McCoy and Wu) of particles being the "mesons" , consisting predom- 
inantly of two quarks confined by a long-range force. Systematic implemen- 
tation of this idea leads to a version of the Bethe-Salpeter equation, which 
yields infinite sequence of meson masses. The Bethe-Salpeter spectrum be- 
comes exact in the limit when the magnetic field is small, and we develop 
the corresponding weak-coupling expansions of the meson masses. The 
Bethe-Salpeter equation ignores the contributions from the multi-quark 
components of the meson's states, but we discuss how it can be improved 
by treating these components perturbatively, and in particular by incor- 
porating the radiative corrections to the quark mass and the coupling pa- 
rameter (the "string tension"). The approach fails to properly treat the 
mesons above the stability threshold, where they are expected to become 
resonance states, but it is shown to yield a very good approximation for 
the masses of all stable particles, at all real values of the IFT parameters in 
the low-temperature regime. We briefly discuss how the Bethe-Salpeter ap- 
proximation can be used to address the case of complex parameters, which 
was the main motivation of this work. 
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1 Introduction 



The Ising Field Theory (IFT) is the quantum field theory in two dimensions which 
emerges in the scaling limit of the Ising model (or any other system from the same 
universality class) near its critical point. The RG fixed point associated with the 
Ising criticality is the c = 1/2 Minimal Conformal Field Theory (see e.g. [1,2]) 
which has two nontrivial relevant operators. Therefore the IFT is understood as 
the Minimal CFT perturbed by the two relevant operators, as described by the 
Euclidean Action 

.Aift = Ac=x/2 CFT + T J e(x)d 2 x + h j a(x)d 2 x , (1.1) 

where the fields e{x) ("energy density") and cr(x) ("spin density") have conformal 
dimensions (5,5) and respectively^- The parameters r and h are scaled 

deviations of the temperature T and the external magnetic field H of the Ising model 
from the critical values, i.e. r ~ T c — T and h ~ H in the limit T c — T — > 0, H — > 0. In 
what follows we refer to h as the magnetic field. It is assumed here that the definition 
of e(x) is such that the domain r > corresponds to the low-temperature (ordered) 
phase of the model. The parameters r and h are dimensionfull, r ~ [length] -1 
and h ~ [length] - ~ . Therefore the physics of IFT essentially depends on a single 
dimensionless scaling parameter, which we define asd 



■n = —r = 2 ^—^r > ( L2 ) 

|n 1 15 I 1 15 

where we have also introduced the notation m = 2ttt which is extensively used 
below. 

Besides its value as the model of near-critical statistical mechanics, the IFT 
appears to be rich a model of two-dimensional particle theory. It was argued in 
the pioneering work of McCoy and Wu [3] that as one changes n from — 00 to 
+00, the spectrum of stable particles undergoes evolution from a single particle to 
an infinite tower of "mesons" formed by weekly confined "quarks". This scenario 
was subsequently confirmed by exact [5] and numerical [4, 6] results, where some 
quantitative details were added. 

1 Here we assume conventional normalization of the fields, such that 

( e(x) e(0) } —> \x\~ 2 and (a(x)a(O) } — > |a;p 3 as x — > . 
This condition fixes the normalizations of the couplings r and h. 

2 This is the same definition of r\ as was used in [4,22] ; note that it differs from definition 
in Refs. [16,29,30]. 
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The IFT admits an exact solution in few special cases. First, at zero magnetic 
field, i.e. at 77 = ±00, the IFT is well-known to be the theory of free Majorana 
fermions (see e.g. [1,7]) of mass |m|. Second, the IFT is shown to be integrable 
at m = and h 7^ (i.e. at 7/ = 0) [5]. In the last case, the associated particle 
theory involves eight stable particles whose scatterings are described by certain 
factorizable S-matrix. Finally, if one admits complex values of the parameters in 
(jl.ip . yet another integrable theory is attributed to the vicinity of the point (s) where 
h takes special pure-imaginary value h = dbi (0.1893...) m~% . This point is the Yang- 
Lee edge singularity [8]. In the vicinity of this point one of the particles of the 
theory becomes light, and its dynamics is described by the integrable Yang-Lee field 
theory [9-11]. At present there are no indication that IFT (jl.ip might be integrable 
at any other values of the parameters, real or complex. 

This paper is our progress report on (a part of) the program "Ising spectroscopy" . 
The program is about the detailed study of the mass spectrum of the IFT at all values 
of the scaling parameter (with the emphasis on the analytic properties of the masses 
as the functions of 77), using perturbative expansions around the integrable points, 
combined with numerical methods. There are several motivations for investing effort 
in this program. For one, the IFT appears in some ways the most basic unitary 
quantum field theory in two dimensions. The Ising fixed point is the lowest (in 
terms of the central charge) nontrivial unitary CFT in 2D [2], and the theory (jl.ip 
corresponds to generic RG flow which originates at this point. Also, the IFT is a 
playground for studying many phenomena common in QFT, like quark confinement, 
resonance states, first order phase transition and associated "false" vacuum state, 
etc. 

In this paper we discuss the mass spectrum of the IFT at real h and m, in the 
low-temperature (low-T) domain m > 0, i.e. at real 77 > 0. Let us remark here that 
at h 7^ distinction between the low-T and high-T (7/ < 0) regimes is fuzzy, and 
to a large extent artificial. Physical quantities like the particle masses pass through 
the point 77 = continuously (in fact, analytically). However, although continuously 
related, the two regimes call for an emphasis on different phenomena, and often 
require different theoretical treatment. Therefore the distinction is conventional in 
the literature on the subject, and we will usually hold to it here and in forthcoming 
reports on our program. 

According to McCoy- Wu scenario, at rj < 00 the free fermions of the h = theory 
(the "quarks") become confined in "mesons". The spectrum of the mesons is dense 
at large 77. As 77 decreases, the heavier mesons become unstable against decay into 
the lighter ones, and successively disappear from the spectrum of stable particles, 
typically becoming resonances. As 77 approaches zero, there are only three stable 
mesons (and many resonances) left (see [4]), and at 77 = some of the resonances 
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Figure 1: The mass spectrum M n (rf) of IFT particles at positive (and some negative) rj. 
The solid lines represent numerical data obtained using the TFFSA. The dotted line shows 
the stability threshold 2M±; after crossing this line the particles become unstable, and their 
masses (not shown) develop imaginary parts. The exception is the point rj = 0, where there 
are eight stable particles, whose masses are indicated by bullets • . The dashed lines show 
the masses obtained from the Bethe-Salpeter equation (|4.8p , with the renormalized string 
tension / defined as in Eq. (|8.3p (see Section 8 for explanations). The circles o indicate 
positions of the higher thresholds Mi + M2, -Mi + M3, etc, at 77 = 0. 



regain stability (the phenomenon related to integrability of the IFT at 77 = 0), so 
there are exactly eight stable particles present[l. Very accurate numerical data for 
the particle masses M n at all 77 can be obtained using the "Truncated Free Fermion 
Space Approach" (TFFSA)B; the results for positive r\ are shown in Fig.l. The form 

3 The process of depletion of the spectrum continues when 77 crosses to the negative values, until 
below a certain point (77 < 772 ~ —2.09, see [4]) there is only one stable particle left. At r/ — > —00 all 
interactions disappear and this particle becomes a free one, in accord with the IFT at zero h being 
a theory of free fermions. We have much more to say about detailed structure of the spectrum in 
this domain, but will deffer this discussion to another paper. 

4 The TFFSA is a modification of the "Truncated Conformal Space Approach" (TCSA) of Ref. 
[12], particularly suited for the IFT, see Ref. [4] for details. 



3 



of the mass spectrum at sufficiently large r? is consistent with the McCoy- Wu idea 
of the particles being composed of two weakly confined quarks. This idea allows 
one to develop the weak-coupling (large- rj) expansions of the masses M n {rj). The 
leading term was first obtained in the original paper [3], and several further terms 
added in Ref. [4]. On the other hand, for highly excited mesons the semiclassical 
approximation applies, and the corresponding expansion was recently derived by 
Rutkevich [13] El- 

In this work we develop a more systematic approach based on the two-quark 
approximation. In Sections 3 and 4 we derive the Bethe-Salpeter (BS) equation, i.e. 
the relativistic two-quark bound state equation for the IFT mesons. It is not exactly 
a textbook exercise because no covariant perturbation theory for the IFT is available, 
and one has to do with the ordinary non-covariant approach. In Section 4 we argue 
that the Lorentz covariance is restored by adding certain corrections to the quark 
self-energy, which technically come from the multi-quark sectors. The resulting BS 
equation (Eq. (|4.8|) ) is very similar to the bound-state equation previously obtained 
by 't Hooft for the two-dimensional model of multicolor QCD [14], which is not 
surprising because the physics is similar. The two-quark approximation can be 
improved by taking the multi-quark sectors into account perturbatively. In Section 
6 we discuss the role of the multi-quark corrections, in particular the effects of the 
renormalizations of basic parameters - the quark mass and "string tension" - due 
to the multi-quark effects. In Sections 7 and 8 we find the spectrum of the BS 
equation numerically, and make a tentative comparison with the TFFSA data. If 
the renormalization of the string tension is taken into account, the BS spectrum 
reproduces the actual masses rather accurately at all positive rj, as is shown in 
Fig.l. The actual masses M n (rj) exhibit a subtle behavior^ when they approach the 

5 We have worked out the semiclassical expansion independently, before [13] appeared. Our 
approach is described in the Section 5 below. We would like to stress that beyond the leading order 
our result (|5.25|l differs from the corresponding term found in [13]. The discrepancy can be traced 
to the corrections to the quark self-energy, which we believe are properly taken into account in the 
Bethe-Salpeter equation (|4.8|l . but are missing in the treatment in [13]. 

6 The near-threshold behavior of the meson masses M n (rf) deserves a separate discussion. Pre- 
liminary analysis (which will be reported elsewhere) shows the following general pattern. The n-th 
meson exists as a stable particle while n exceeds certain threshold value, n > r)„, in which domain 
M n (rj) remains below 2Mi(n). At 77 = rj n the graph of M n (rj) touches the curve 2M\(rf), so that 
2M\{rf) — M n (rf) ~ (r/ n — rj) 2 as 77 — > rj n . Correspondingly, the analytic continuation of AI n (rj) to 77 
immediately below r/ n remains real, although no stable particle with this mass exists at 77 < rj n (the 
situation known as the "virtual level", see e.g. [15], §128). By further decreasing the value of rj one 
reaches singularity of the function M n (r\) (typically a square-root branching point) at rj = rj n < rj n , 
and below f\ n the mass M n (rj) develops an imaginary part, so that the n-th meson reemerges as a 
resonance state. This pattern seems to hold for all M n with n > 3. Our preliminary estimates of 
the thresholds are 772 ~ —2.06, 773 ~ —0.136, 774 ~ 1.0, while 773 m —0.515, and the higher fj n are 
very close to 77„. 
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"stability threshold" 2M\ (it is visible in Fig.l for the mass M3, and less prominently 
for M4 and M5), the feature certainly driven by the multi-quark components of the 
meson states. The Bethe-Salpeter approximation does not capture these subtleties. 
Also, the BS masses remain real above the stability threshold, while the actual 
IFT particles in this domain generally become unstable, and their masses develop 
imaginary parts. Although in this paper we do not discuss the resonance part of the 
IFT spectrum in any detail (we plan to address this important subject separately; 
some results on the IFT resonances can be found Ref's [6, 16]), let us note that 
the pattern shown in Fig.l suggests that when 77 goes to zero some of the the IFT 
resonances turn into the particles of the corresponding integrable theory, while others 
become weakly coupled bound states of those, with the binding energy ~ if. For 
instance, 

M 5 (r/) = Mi (r ? )+M 2 (r ? ) + 0{ V 2 ) as rj -> . (1.3) 

This is one of many interesting phenomena which the IFT spectrum appears to 
exhibit; it provides illustration of its remarkable richness, and gives us additional 
motivation for its systematic study. 



2 Ising quarks and Ising mesons 

At h = 0, and in the low-temperature phase m > 0, the spin-reversal symmetry of 
the IFT is spontaneously broken, and the theory has two degenerate ground states 
differing by the sign of the spontaneous magnetization 

(a(x))± = ±a , (2.1) 

where a is a constant whose exact value is [17] 

a = \m\*s, s = 2 1 / 12 e-t^'^ 1 ) = 1.35783834170660... (2.2) 

At the same time, as was mentioned in the Introduction, at h = the IFT (11. ip 
reduces to the theory of free Majorana fermions of the mass m. These free fermions 
are identified with the domain walls separating the spatial domains of positive and 
negative magnetization. Such interpretation makes it obvious that adding the in- 
teraction term h J a(x)d 2 x with small h (h « m~a ) to the free theory generates 
confining force between the domain walls, with the "string tension" 2ah. Due to 
this effect, the stable particles appear as the "mesons", i.e. the bound states of two 
"quarks" 

7 The gauge group associated with this confining interaction is Z2 (the Ising model with magnetic 
field is known to be dual to the Z2 gauge theory, see e.g. [1]), hence the quarks and the antiquarks 
are identical. 



5 



The above picture alone allows one to give some quantitative description of 
the masses M n , n = 1,2,3,... of the mesons in the limit rj — ► +00. Indeed, the 
Hamiltonian of two particles interacting via the confining force is 

H = uj(pi) + uj(p 2 ) + 2ah \xi - x 2 \ , (2.3) 

where xi,x 2 and pi,p 2 are the coordinates and the momenta of the quarks, and 

w(p) = ^m 2 +p 2 . (2.4) 

For the lower part of the meson spectrum {M n with n fixed) the non-relativistic 
approximation uj(p) ~ m+p 2 /2m applies, and the problem reduces to finding energy 
levels of the Hamiltonian p 2 /m + 2ah\x\. Since the quarks are fermions, only odd 
levels are relevant, and one finds 

M n -2m -► (2ah)* z n as h -► (2.5) 

where z n ,n = 1,2, ... are consecutive zeros of the function Au(— z). This result was 
first obtained by McCoy and Wu [3] in their analysis of the spin-spin correlation 
function. For higher mesons (M n with n ~ m 2 /ah) the semiclassical approximation 
is more suitable. Classical trajectories of the system (??) are periodically repeated 
cycles in which, during the time between two collisions, the quarks move under 
constant acceleration directed towards each other. It is convenient to parametrize the 
time t as (t — to) = R sinh/3, where to is suitable reference time, and R = m/ (2ah). 
The parameter f3 has a simple meaning: +/3 and —f3 are the rapidity of the two 
quarks in the center-of-mass frame. In this frame, the classical trajectory within 
each cycle is given by the equations (up to permutation of the quarks) 

p = ^ 1 P 2 ^ = _ m gjnh p ^ x = x\ — x 2 = 2R[ cosh d — cosh (3) (2.6) 

for — 1} </?<#, where $ is positive parameter characterizing the classical trajectory. 
At P = ±$ the separation between the quarks vanishes, hence the mass (i.e. the 
center-of-mass energy (??)) associated with the trajectory is M = 2m coshd. The 
reduced action per cycle is 

[ pdx = 2mR f sinh 2 (3d(3= sinh2 ^~ 2 ^ . (2 . 7 ) 

Jcycle J— d ^ 

here and below we use the notation A for the dimensionless ratio 

A = ??*. < 2 -8> 

m A 
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If one first treats the quarks as distinguishable particles, the full period consists 
of two cycles, so that the reduced action per period is twice the value of (12. 7j) . so 
that the Bohr-Sommerfeld quantization condition reads 2 f cycle pdx = 2tt(N + 1/2) , 
where the integers N must be taken to be odd, N = 2n — 1, to reclaim the fermionic 
nature of the quarks. The quantization condition then takes the form 

sinh 2i9 n — 2-Q n , , s , 

\ = 2vr(n- 1/4) , n= 1,2,3,... (2.9) 

A 

leading to the WKB mass spectrum of the mesons, 

M n = 2m coshtf n . (2.10) 



3 The two-quark approximation 

The idea of the Ising particles being predominantly two-quark composites allows one 
to develop somewhat a more systematic theory of mesons, in particular to develop 
a weak-coupling expansions for their masses M n . One starts with the free-fermion 
theory at h = 0, and introduces the creation and annihilation operators aj and a p 
(where p denotes the spatial momentum of the particle), subject to the canonical 
anticommutators {a p ,aj} = 2ir5(p — q). Below we will use the notation 

\Pl,P2,-~ ,Pn) = a pi a P2 ' - '4n I ) ( 3 - 1 ) 

for the n-particle states. 

The h = theory is described by the free Hamiltonian 

u(p) a£a p JL , (3.2) 

-oo 27T 

where u)(p) is the relativistic energy (|2.4|) , and Eq is the ground-state energy. The 
interaction is generated by adding the term associated with the last term in the 
action (jl.ip . 

/oo 
<t(x) dx , (3.3) 
-oo 

where x is the spatial coordinate, and <r(x) = cr(x, t) |t=o ■ The operator <r(x) can not 
be expressed through the fermions in a local way (see e.g. [1]), but its matrix elements 
between the states with any number of the fermions are known explicitly [18]. Here 
we will only need the 2 — > 2 matrix element 

(pi,p 2 | o-(O) \qi,qz) =4a g(p 1 ,p 2 \qi,q2) , (3.4) 
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where 



G(Pi,P2\qi,Q2 



(1/4) 



\/ w (Pl)w(P2)w(«l)w(?2) 



u(pi) + w(gi) w(p 2 ) + u(q 2 



+ 



Pi - 9i 
Pi -P2 



P2 ~ q2 
qi - q2 



Pi -92 P2-51 + ^(p 2 ) w(?i) + w(g 2 ) 

One looks for the eigenstates of the Hamiltonian (|3.3p of the form 



$(2) \ + I $(4) \ + I ^(6) v + 



where 



*(2) 



dpi 
2vr 



dp_2 
2vr 



*(pi , Pa) I Pi,P2! 



(3.5) 



(3.6) 



(3.7) 



is the two-quark component, and further terms in (|3.6j) represent contributions from 
the multi-quark (i.e. the four-quark, six-quark, etc) sectors. The two-quark approx- 
imation is developed under the assumption that (|3.7j) is the dominating component 
of the meson state. It is certainly valid at sufficiently weak coupling, when the pa- 
rameter A is small, but we will see that it provides meaningful description of the 
spectrum even at very large values of this parameter. Anyway, the two-quark ap- 
proximation can be improved by taking the multiquark components in f)3.6[) into 
account perturbatively (see Sect. 6). 

The two-quark component is completely characterized by the "meson wave- 
function" VP(pi,2>2) m (13. 7|) . By the definition (|3.7p the wave function is antisymmet- 
ric, Vl/(pi,p2) = — ^(P2,Pi)- If the multi-quark components in f|3.6|) are neglected, 
the eigenvalue problem (H — E) \ Sifp) = reduces to the integral equation 

[e(pi) + e(p 2 ) - AE]y{ Pl ,p 2 ) = 

poo 

f f 2TT5{p l +p 2 -q 1 -q 2 ) g(p llP2 \ qi , q2 )*( qi , q2 )^^ (3. 8 ) 

J — oo 



2. 



where at this point e(p) stands for the free quark energy e{p) = u>(p) = \J m 2 + p 
we introduce this new notation here in preparation to the discussion in Sect. 4 below, 
where e(p) will become the "dressed" quark energy. In ()3.8f) AE = E — E vajC is the 
energy above the ground state, and 



fo = 2ah = m 2 A 



(3.9) 



is the "string tension". In fact we will refer to the parameter /o, Eq. (|3.9p . as the 
"bare" string tension, as opposed to "dressed" , or "effective" string tension / which 
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will replace fo when radiative corrections (originating from the multi-quark com- 
ponents in (|3.6j) ) are taken into account, as will be discussed in greater details in 
Sect. 6. The kernel Q in the right-hand side is the matrix element (|3,5p . The kernel is 
singular, and the r.h.s. of (|3.8p involves the principal value of the singular integral. 
For the meson state with the momentum P one takes the wave function of the form 

*(pi,p 2 ) = (27r)*(pi+p2-P)*p( P i-P/2) (3.10) 

with antisymmetric *$>p(p) = —^>p(—p). Then the equation (|3.8p takes the form 

[e(P/2-p)+e{P/2 + p)-AE}* P (p) = f f Q ° G P (p\q)* P (q) ^ , (3.11) 
where the kernel Gp(p\q) is the function (|3.5p evaluated at p\ + p2 = qi + <?2 = P, 



f + f '/ ) ■ ( ; >-J- ) ) 



The meson wave- function ^?p(p) is assumed to be normalizable H, 



and the equation (|3.1ip is understood as the eigenvalue problem for the meson energy 
AE. 

Observe that the kernel Gp{p\q) in (|3.1ip has the second-order poles at q = ±p, 
with the residues 1. More precisely, one can check that 

G P (p\q) = -^— 2 - — ^ + G%^\p\q) , (3.14) 
[P - Q) (P + 1) 

where the last term is regular when both p and q take real values. Therefore the prin- 
cipal value integral in the r.h.s. of (|3.11|) has the same effect as the linear potential in 
(|2,3p . provided the separation between the quarks is large enough, m |xi — X2I >> 1. 



The short-range interaction described by the regular term G^p eg \p\q) is what makes 
(|3.1ip different from the simple-minded Eq. (|2.3p . 

The integral equation (13. lip is the result of straightforward application of the 
idea of the meson being the two-quark construct. However, it is well known that in an 



8 More precisely, the states (|3.7[) are assumed to be 5-normalizable in terms of the meson mo- 
mentum P. Since 

(*p I *p,) = 2ttAE5(P-P') ||*p|| 2 

this implies the Eq. (|3.13p . The factor AE is inserted here in anticipation of relativistic normaliza- 
tion of the meson states. 
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(a) 



(b) 



Figure 2: Possible world lines of quarks in a meson, (a) Both quarks propagate forward in 
time, (b) Creation an annihilation of virtual pairs leads to the presence of more then two 
quarks in the intermediate state. 

interacting theory neglecting the multi-quark components of the state (|3.6p generally 
violates Lorentz invariance. This is because in a relativistic theory the quarks must 
be allowed to move occasionally backward in time, the possibility depicted in Fig. 2 
but obviously neglected when all the multi-quark terms in (|3.6p are discarded. As 
the result, the equation (|3.1ip is not Lorentz-invariant, in particular the eigenvalues 
AE are not expected to have correct dependence V M 2 + P 2 on the momentum P. 

To restore the Lorentz invariance, it is necessary to retain certain contributions 
from the multi-quark sectors in (|3.6|) , the ones associated with the backward- in-time 
propagations of the type depicted in the Fig.2b. Adding such contributions would 
result in the Lorentz invariant version of the two-quark approximation. Analogous 
problem in perturbative quantum field theory is solved by the Bethe-Salpeter equa- 
tion (see e.g. [19]). Unfortunately, in the IFT manifestly covariant perturbation 
theory in h is not available. For that reason, in the next Section we discuss the 
multi-quark contributions associated with the backward-in-time propagation of the 
quarks. These contributions are described as corrections to the dispersion low of 
the quarks inside the meson. We will argue that these corrections disappear in the 
limit P — > oo, and therefore the Lorentz-invariant Bethe-Salpeter equation is essen- 
tially the equation (|3.1ip in the infinite-momentum frame. Similar arguments were 
previously put forward in the context of the 'tHooft model in [20] 
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(a) 



(b) 



(c) 



Figure 3: Diagrams of the perturbation theory in h. (a) The "ladder" diagrams. The blobs 
represent the matrix elements (|3.4I) . These diagrams are summed up by the Eq. (|3.11[) . 
(b) The diagram with four intermediate quarks, (c) Example of a disconnected part of the 
diagram in (b). 

4 The Bethe-Salpeter equation 



As usual, the integral equation (|3.1ip can be understood as a device for summing 
up certain class of diagrams of the perturbation theory in h. Here and below we 
have in mind usual quantum-mechanical perturbation theory (see e.g. [15]), not the 
covariant Feynman perturbation theory □, so that all the intermediate lines in the 
diagrams below are on-shell. The equation (|3.1ip takes into account the "ladder" 
diagrams shown in Fig. 3a, where the blobs represent the matrix elements (I3.4p . with 
only two quarks involved in the intermediate states. But the Eq. (|3.1ip neglects all 
the diagrams with four, six, or more intermediate quarks, like the one shown in the 
Fig. 3b. The blobs in the Fig. 3b stand for the matrix elements between the two- 
quark and four-quark states, and the integration over the intermediate momenta 
ki,k2,k^, &4 (subject to the constraint k% + k 2 + ^3 + k^ = pi + P2 = Qi + Q2) is 
implied. The contribution of the diagram in Fig. 3b contains disconnected parts, 
in which only one of the particles is affected by the interaction, as is depicted by 
the diagram in Fig. 3c. In the disconnected parts like this the momentum conserves 
separately within each disconnected component, i.e. such diagrams contain two 
momentum 5-functions. In particular, the contribution of the diagram in Fig. 3c to 
the energy can be written as 



disconnected 
diagram in Fig. 3c 



(2vr) 2 5( Pl - qi )5(p 2 - q 2 ) 4 3) (P2) 



(4.1) 



9 As was already said, manifestly covariant perturbation theory in h is not yet developed. See 
however calculations in Ref. [21] which seem to be suggestive in this respect. 
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l 2k+l 





(a) 



(b) 



(c) 



Figure 4: The quark self-energy: (a) The simplest self-energy diagram of the order ~ h 2 ; 
its contribution Y$p(p) appears in the Eq. (|4.1[) . (b), (c) Examples of more complicated 
self-energy diagrams. 



where E^Cp) i s the self-energy term from the diagram in Fig. 4a. Clearly, (|4,ip 
and similar disconnected contributions involving the self-energy diagrams in Fig. 4 
can be accounted for by replacing the factors uj(p) = \J m? + p 2 in the left-hand 
side of the Eq. (l3.8p by the dressed quark energy, uj{p) — > e(p) = ui(p) + E(p), 
where the self-energy correction in principle can be determined order by order 
in the perturbation theory in h. On general grounds, one expects to have two kinds 
of perturbative contributions to the self-energy. As usual, there are contributions 
associated with the one-particle irreducible diagrams. These corrections do not 
modify the relativistic form (|2.4p of the momentum dependence of the quark energy, 
while making corrections to the quark mass m. At this stage we shall ignore the 
radiative corrections to the quark massif, postponing the discussion to Section 
6. Here we are more interested in the contributions from the one-particle reducible 
diagrams. Simplest example of the one-particle reducible diagram is shown in Fig. 5b, 
which can be interpreted as the quark traveling backward in time in between two 
interaction events. Note that unlike the one-particle reducible diagram with the 



The leading correction ~ h to the quark was determined in Ref. [22] (the result is quoted in 



Eqs. (|6.2p . (|6.3|l below). It may be useful to note that the calculations involve subtracting the one- 
particle reducible parts (Eqs. (5. 10), (5. 13), and Fig. 3 of Ref. [22]), which are exactly the diagrams 
Fig. 5a and Fig. 5b here. The leading correction in (14. 4[) below comes from adding the diagram in 
Fig. 5b back. 
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Figure 5: One-particle reducible self-energy diagrams. Parts (a) and (b) show two types 
of the order h 2 diagrams, with one and three quarks in the intermediate state, respectively. 
Part (c) shows an example of the higher-order diagrams which contribute to s(p). 

forward time propagation in Fig. 5a, which is already accounted for in the sum of 
the ladder diagrams in Fig. 3a, the diagram in Fig. 5b has to be taken into account 
as separate contribution to e{p). The difference between the diagrams Fig. 5a and 
Fig. 5b is in the time order of the interaction events; since this distinction may depend 
on the Lorentz frame, the contribution of the diagram in Fig. 5b is not expected to 
be Lorentz covariant. Similar non-covariant corrections to the quark self-energy 
arise from the one-particle reducible diagrams of higher orders with the backward- 
in-time propagation (e.g. the diagram in Fig. 5c). The quark energy e(p) includes 
all the one-particle reducible diagrams which are not reproduced by iterations of the 
equation (|3. (such are the diagrams in Figs. 5b and 5c, but not in Fig. 5a). 

It is instructive to calculate explicitly the lowest-order correction term in e(p), 
This term comes from the diagram in Fig. 5b, which reads 



Performing the trivial integration, and using explicit form [18] of the matrix element 






(4.3) 



one finds 




(4.4) 
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where A is defined in Eq. (|2.8p . It is possible to show that this correction in (|4.4p is 
exactly what is needed in order to restore the Lorentz covariance of the spectrum of 
the equation (|3.8p , to this order in A. The calculation supporting this statement is 
rather involved; we will present it elsewhere. 

In principle, the higher order terms in (|4.4p can be determined by direct calcu- 
lation of the corresponding diagrams. It is reasonable to expect that incorporating 
these terms in e(p) will restore the Lorentz covariance of the equation (j3.8l) order 
by order in A 2 (analogous statement for the 't Hooft's model has been proven in 
Ref. [20]). In particular, with all contributions to taken into account, the 

P-dependence of all the eigenvalues A.E of (|3.8p must assume the relativistic form 

AE = \JM 2 + P 2 (4.5) 

with M interpreted as the meson mass. Let us stress that at the moment we have 
no proof of this statement beyond the order A 2 . We intend to come back to this 
question in the future. In the remaining part of this paper we proceed under this 
assumption. 

Once the assumption is accepted, there is an obvious way to circumvent calcu- 
lation of S(p) altogether. Note that the leading term ~ A 2 in (|4.4j) decays fast (as 
|p|~ 3 ) when p — » oo. It is possible to verify that the higher-order terms in (14. 4p 
decay yet faster, so that the asymptotic form 

£ ^ = ^ + ^] + (lpj) ( 4 -6) 

is valid to all orders in Ap^i Consider the eigenvalue equation (|3.1ip and take the 
limit P — ► +oo. As is shown in Appendix A, in this limit, and under appropriate 
normalization, the meson wave function ^p(p) remains finite at 2\p\/P < 1, but 
vanishes as 1/P 2 outside this domain. For 2\p\/P < 1 both sides of the equation 
(13. lip decay as 1 jP at large P. Note that this leading term of the P — > oo asymptotic 
is not affected at all by the correction terms in (|4.4p . since S(p) decays faster then 
1/p. Balancing the coefficients in front of 1/P asymptotic in both sides of (|3.1ip 
yields a nontrivial eigenvalue equation for the parameter M 2 in (|4.5p . This equation 
is derived in Appendix A. Its most convenient form is obtained by changing to the 
rapidity variables 

2p = P tanh 9 , 2q = P tanh 6' (4.7) 

11 Here we ignore all corrections to the quark mass determined by the one-particle irreducible 
self energy diagrams. Of course, exact asymptotic form of e(p) is given by the Eq. (|4.6[l , with m 
replaced by the "dressed" quark mass m q . 
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before sending P to infinity, and introducing the notation ^(9) for the limit of ^p{p) 
at P — > oo. It is 



M 2 



m 



4 cosh" 1 9 



/OO Jfl/ 
G^')^) — , (4. 



where 



G{9\9') 



2 cosh(0-0') I sinb6» sinhfl' 



sinh 2 (6> - 6>') 4 cosh 2 cosh 2 0' 



(4.9) 



We will refer to (14.81) as the Bethe-Salpeter (BS) equation. 

Let us make few general remarks on the properties of the BS equation (|4.8p . 
First, the wave function ^>{9) should be regarded as a vector in the Hilbert space 
with the metric 



1*1 



|<f(0)| 2 d9 



4 cosh 2 9 2vr 



(4.10) 



which is just the metric (|3.13p rewritten in terms of the rapidity variable 9, in the 
limit P = oo. Second, the equation is understood as the eigenvalue problem for the 
parameter M 2 , 



where the operator H, defined as 



H^(9) = 4 cosh 2 9 



m 2 (0) - f 



rlf) 

G(9\9>)*(9>)- 



(4.11) 



(4.12) 



is Hermitian with respect to the metric ()4.10p . Its eigenvalues are real and positive. 
We will use the notations M 2 with n = 1,2,3, .. . for the successive eigenvalues of 
([4. 11|) (arranged in the order of increasing magnitude) and ^ n {9) for the correspond- 
ing eigenfunctions. The quantities M n provide certain approximations for the actual 
meson masses, and we reserve the notation M n for those actual masses. 



5 Weak coupling expansions 

If the magnetic field h is small, such that /o << m 2 , the Bethe-Salpeter equation 
(14. 8p reproduces physics of weakly confined quarks, as described in Sect. 2. According 
to this picture, one expects that in the limit A — ► each individual eigenvalue M 2 
approaches 4m 2 from the above, so that the full spectrum {M 2 } of the equation 



15 



(|4.1ip becomes dense in the segment [4m 2 , oo). Here we develop the weak coupling 
expansions of the eigenvalues M 2 . The nature of the expansion depends on what 
part of the spectrum {M 2 } one is interested in. First, one might be interested in the 
small- A expansion of an individual eigenvalue M 2 , with fixed n. Since the kinetic 
and potential energies of the quarks in any given meson state become small at small 
A, we refer to this case as the low energy expansion. The equation (|2.5p is just the 
first therm of this expansion, and in general the low-energy expansion of M 2 turns 
to be in fractional powers of A, namely in the powers of the parameter 



i 



t = A» . (5.1) 

Second, the behavior of the higher levels M 2 with n > i may be of interest; in this 
case one can derive the semiclassical expansion extending the Eq. (12.10p . This ex- 
pansion is in integer powers of A, and the Bohr-Sommerfeld equation (|2.10p appears 
as its leading term, while the higher-order corrections were discussed in [13]. In all 
cases it is useful to remember that the BS equation itself is an approximation, and 
therefore starting from certain order (i 9 in the low-energy and correspondingly A 3 
in the semiclassical expansion) both expansions exceed the accuracy of the BS equa- 
tion (|4.8p . From these orders on, the coefficients of the weak coupling expansions of 
the actual masses M n start to receive contributions associated with the multi-quark 
components of the meson states^! . In this Section we ignore such corrections, but 
we will come back to them in Sect. 6. 

We develop the weak coupling expansions starting from the following approxi- 
mation for the solution of the equation (|4.1ip , 

t,(ow • x.6 f°° eim cosh /3 dp , r n , 

* (0) (*) = sum J ^ s . mHe + p _ i0) smh{e _ p + l0) ■ (5-2) 

where 

M 2 

S(P) = —r tanh p - . (5.3) 
4?7r 

It is easy to show that is normalizable (with respect to the metric (|4.10p ) at 

any value of M 2 . The quality of this approximation is determined by the function 
A(6\M 2 ) defined as (in what follows we usually suppress the argument M 2 ) 

Am 2 sinhfl A(0) = [H - M 2 ] 9 (0) , (5.4) 



12 Here we refers to the multi-quark corrections which cannot be absorbed into the renormaliza- 
tions of the parameters m and fo- Corrections to these parameters appear at lower orders; more 
details are in Section 6. 
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where H is the operator (|4.1ip . By straightforward transformations (which we sketch 
in the Appendix B) this function can be written as the integral 

/oo 
-oc 

where 



cosh 2 fl / (9 + 13) (9-13) \ 1 g_ l 

sinhfl \swh.(d + (3) swh.(d-(3)) 4 cosh/^J ' { ' ' 

Note that the terms ~ A in the integrand in (|5.5p are sorted out according to their 
symmetry with respect to [3 — ► — /3. 

Next, we develop a small -A expansion of the integral (|5.5p . Let us note that 
the integrand in (|5,5p is independent of # except for the last term B(9\(3) in the 
brackets there. If the 9 dependence could be ignored, it would be possible to turn 
the whole integral (|5.5p to zero by tuning the parameter M 2 to special values, the 
zeros of the function A(M 2 ) = A(9\M 2 ); these special values would then determine 
the eigenvalues M 2 . It turns out that the term B(9\(3) in (|5.5p does not contribute 
to the expansion all the way up to rather high order - up to t 8 of the low-energy 
expansion and up to A 2 in the semiclassical expansion - and therefore evaluating 
first few terms of the weak coupling expansion from (|5.5p is very easy. At higher 
orders the ^-dependence shows up and hence no choice of M 2 can turn the function 
(15. 5p to identical zero. This is well expected - after all (|5.2p is an approximate, 
not exact solution of the Eq. (j4.8j) . Therefore finding the terms of the order A 3 and 
higher requires additional ingredients. Let us show how to develop systematic weak 
coupling expansion of the eigenvalues M 2 to any orderf* 3 !. 

We start with the trivial observation that any eigenvalue M 2 satisfies the equa- 
tion 

C n (M 2 ) = , (5.7) 

where 

C n (M 2 ) = [H - M 2 ] *(°)) . (5.8) 

Here and below the brackets ( , ) denote the scalar product associated with the 
metric (|4.10p . We also assume the eigenvectors to be orthonormalized, (*$> n ,*$>k) = 

13 In the case under consideration the higher order terms have somewhat limited value: as we 
explain in Section 6 such terms compete with the multi-quark corrections which are generally 
beyond control of the BS equation (14. 8[) . Here we still do it as a technical exercise, since the 
method can be applied in other models, notably in the 't Hooft's model of multicolor QCD [14], 
where very similar Bethe-Salpeter equation (shown as the Eq's (| A.6I) . (I A. 8|) in Appendix A) is exact. 



M 2 
Am 2 cosh (3 



+ 



iX sinh (3 
8~ cosh 2 13 



+ \B(9\P) 



ei s ^df3, 



(5.5) 
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5 n k- Note that 



COSH 6J 2"7T 



At small A the function fy(°\9) provides good approximation for ^ n {6), and more- 
over one can find corrections by iterating the equations 



*n(0) 



(*<°>,*n) 



2\ 



(5.10) 

M 2 =M2 



Then the solution has to be plugged back into (|5.8p . providing further corrections 
to the eigenvalues through the Eq, (|5.7|) , 

We illustrate this technique by explicit calculations of the first few orders, both 
in the low-energy and semiclassical expansions. 

5.1 Low energy expansion 

Assume that the parameter M 2 in (|5.5f) is close to 4m 2 , so that M 2 — 4m 2 < t 2 . 
Note that any given eigenvalue M 2 enters this domain at sufficiently small A. When 
this is the case, the integral (|5.5p is dominated by contribution from ~ t. It 
is convenient then to change to a new variable u = —/3/t . We also write M 2 = 
4m 2 (1 + z t 2 + 0(i 4 )) with some coefficient z to be determined below. In the leading 
order in t the Eq. (|5.5j) reduces to the integral defining the Airy function, 



oo 

* „3 



A(0) = t / ea" ~ lzu du + 0{t 6 ) . (5.11) 



Hence, it is sufficient to set the parameter — z equal to any zero of the Airy function, 

Ai(-z)=0, (5.12) 

to make (|5.1ip vanish to the leading order in t. This reproduces the equation (|2.5p . 
In what follows we use the notation z for a generic solution of the equation (|5.12p . 
More generally, we look for the solutions of (|5.7p in the form 

M 2 °° 

k=3 

The integral (|5.5p can be evaluated order by order in t. Let us write its expansion 
as 

s = D 3 (9) t 3 + D 4 (6) t 4 + D 5 (0) t 5 + • • • . (5.14) 

7T Al (— Z) 
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A straightforward calculation 14 ! yields 

D 3 (9) = -e 3 , D i (9) = ^-e i , D 5 (9) = -e 5 , etc (5.15) 

5 

It turns out that the coefficients D^iO) have no 9- dependence all the way up to 
D$(9), so that they can be turned to zero by suitable choice of €3, €4, ■ ■ ■ , e% , e.g. 
£3 = 0, £4 = z 2 /5, etc. It also turns out that the coefficients £5 and €7 obtained in 
this manner vanish (but higher order odd terms are present in (j5.13fl . see Eq. (j5.19j) 
below). Thus we find 

M 2 2 z 2 4 /3z 3 57 \ fi /23z i 1543z\ „ 

l + zt 2 + — t 4 - — + — W 6 + -_ + __— t 8 + ... (5.16) 



4m 2 5 V 175 280 J V 7875 12600 / 

The first four terms here are identical to the results previously reported in Ref. [4] El, 
but the fifth (~ t 8 ) term is new. Starting with the term ~ i 9 the ^-dependent term 
B{9\j3) in (|5.5h becomes relevant, and the coefficients Dk(9) with k > 9 generally 
have nontrivial dependence on 9. For example 

£> 9 (0) = !! B (9) + ^B 2 {9) + 2B 4 (9)-e 9 , (5.17) 

o o 

where the functions B2 n {9) are the coefficients in the expansion 

B(9\f3) = B (9) + B 2 {9) (3 2 + B 4 (0) /? 4 + • . . (5.18) 

At this order and beyond we have turn to the equation (|5.7|) . and determine the 
coefficients eg, eio, etc from the condition that the integral (|5.9p vanishes order by 
order in t. We do not enter the details of this calculation here, only remark that 
the integral (|5.9p receives the main contribution from the domain 9 ~ t, therefore 
its expansion in t involves the expansions of the coefficients Dg(9), Diq(9), ... in 
the powers of 9 . Let us quote few more terms of the expansion (|5.13p , 

13 1894 z 5 23983 z 2 3313 z 

eq = , ein = , eii = . (5.19) 

y 1120 vr ' 3031875 242550 10080 vr v 1 



'The calculation involves the integrals 

(iu) k e~z u ~ lzu du = 7rAi'(— z) 1^. 



The factors 7^ in the r.h.s. of this equation are easily computed one by one using the recursion 

h+2 = -zh + klk-i with J = and h = 1 , 

where the last two equations follow from Eqs. (|5,12[) (|5.15|l . 

15 Our calculations in Ref. [4] were performed using the Eq. (|3.11[) . with e(p) = lo(p), in the center- 
of-mass frame P — 0. Corrections to e(p) contribute to M„ starting only from the order t s ; that is 
why all terms quoted in [4] are correct. 



19 



Quick inspection of (|5.16p and (|5.19p (as well as the higher order terms) reveals 
a general structure. The coefficients are polynomials in z, of growing degree. The 
even coefficients have the form €21 = ai z l + h z l ~ 3 + • ■ • , while for the odd terms 
we have 621+1 ~ z l ~ 4 . The parameter z becomes large for higher eigenvalues (the 
Airy zeros grow as z n ~ n»). If t « 1 but z n t 2 ~ 1 or greater, it is useful to sum 
up the leading terms ~ {nit) 21 . It is not hard to see that this is exactly what the 
Bohr-Sommerfeld equation (|2.9p . (|2.I0p accomplishes. Then one can try to collect 
the subleading terms ~ i 6 (n3 t) 21 , ~ i 9 (ns t) 21 , etc. This is done by evaluating 
the higher order terms of the semiclassical expansion, as we discuss next. 



5.2 Semiclassical expansion 

If A << 1 and M 2 is not too close to Am 2 the integral (|5.5p can be evaluated by 
the stationary phase method. There are two saddle points at (3 = ± t?, where $ is a 
positive parameter defined by the equation 

M 2 = Am 2 cosh 2 •& . (5.20) 

Correspondingly, the integral can be written as 

A(0) = A+(0) + A_(0) (5.21) 

where the two terms represents the contribution of the two saddle point (3 = ± i9, 
respectively. Obviously, A_(#) = A+(#), so we can concentrate attention on one 
of this terms, say A + (9). In the leading approximation only the first term in the 
square brackets (I5.5P contributes, and it has to be evaluated at the saddle point 
(3 = Thus in the leading order A(#) is independent of 6, namely 



, = = cos — . (5.22) 



where 



gW S ^)U = Si " h2 "- 2tf (5-23) 

(5(t?) is not the same function as S("d) because -Q enters also through ([5720])). Van- 
ishing of (|5.22p leads exactly to the quantization condition (|2.9p . 

Corrections to the leading WKB formula can be obtained by developing the 
standard loop expansions of the integral (|5.5|) around the saddle points. In the one- 
loop order the ^-dependent term B(6\f3) still is not too important (its effect is in 
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bringing in an overall factor 1 + XB(9\'&) in A(0)), and A(9) still can be turned to 
identical zero by appropriate adjustment of Explicit calculation yields 




= e 4 e A 



^W(l + a5 1 ( 1 9)+0(A 2 ) 



(5.24) 



\/27rAsinh2t? 1 + \B(6\0) 



where 



i r 5 i ii ii 



— sinh 2 $ 



(5.25) 



sinh 2# 24 sinh 2 ■& 4 cosh 2 ■& 12 6 



This leads to the one-loop correction to the quantization condition (|2.9p 



sinh2i? n - 2tf n = 27r A (n — 1/4) — A 2 Si(tf n ) - <3(A 3 ) . 



(5.26) 



The eigenvalues M 2 



are still related to through (|5.20p . i.e. 



M 2 = 4m 2 cosh 



n 



(5.27) 



Let us remark here on the relation to the results of Ref. [13], where the semiclas- 
sical expansion to the same order ~ A 2 is developed on the basis of the two-quark 
approximation in the center of mass frame. Our result (15. 26ft is similar but different 
from that of [13]. The calculation in [13] ignores the correction ~ A 2 to the quark 
self energy, which is responsible for the Lorentz covariance of our approach. It is 
easy to check that when the self-energy correction (I4.4p is added, the result of [13] 
becomes identical to (|5.26|) . (|5.27|) above. 

The higher-loop corrections generally leave the Eq. (|5.27p unchanged while adding 
the higher-order terms 



to the r.h.s. of the Eq. (|5.26p . Calculation of the coefficients &2 ($), 5*3 (#),•• • requires 
(rather straightforward) evaluation of the higher-loop contributions to A + (#) and 
A_(0) in (|5,2ip . At two loops (and of course at all higher orders as well) A(0) 
acquires an essential dependence on 8, which makes it impossible to turn it to 
identical zero by an adjustment of $. Just like in our previous analysis of the 
low-energy expansion, one finds the higher-loop terms in (I5.28P from the equation 
(|5,7p . Since A(0) remains slow- varying function at small A, the integral in (|5,9p 
is dominated by the stationary-phase points 9 = ±i? of the function ^f n (9), and 
can be evaluated order by order in the associated loop expansion. Thus, at the 
two-loop level in the mass spectrum calculation, the leading-order stationary-point 



A 3 S 2 ($ n ) - A 4 SM 



(5.28) 
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calculation of (|5.9p is appropriate, therefore the coefficient S2O?) i n f|5.28j) is found 
from the equation 

A (tw °- loop) (0)| e=)? = . (5.29) 

Calculations of yet higher coefficients 6*3(1?), £4(1?), etc involves also iterations of 
(|5,10|) . along with the higher- loop evaluation of the integral (|5.9p . We did not per- 
form these straightforward calculations explicitly because, as was mentioned above, 
all terms beyond A 2 in (|5.26p exceed the accuracy of the BS equation itself - starting 
from the term A 3 the expansion of the actual masses M n receives nontrivial contri- 
butions from the multi-quark sectors, and without good control of these corrections 
the semiclassical analysis at the order A 3 and beyond does not seem to make much 
sense. 

To conclude this section let us remark that the weak coupling expansions (both 
the low-energy and semiclassical) of the meson masses are only asymptotic ones. For 
the expansions of the BS eigenvalues M n this fact can be deduced from their ana- 
lytic structure as the functions of A: there are infinitely many complex singularities 
accumulating at the point A = 0, as we explain in Sect. 9. It is possible to argue that 
the point A = is an essential singularity of the actual masses M n as well, but this 
discussion goes beyond the scope of the present paper. 

6 Multi-quark corrections 

As was already mentioned in the Section 4, the Bethe-Salpeter equation ignores 
substantial part of the multi-quark contributions to the meson masses. In general, 
explicit computation of such contributions is a rather difficult task which we do 
not attempt to perform here. Instead, we briefly discuss the role of the multi- 
quark corrections, in particular the way they modify the weak-coupling expansions 
described in the previous Section. 

The multiquark corrections are represented by diagrams with more then two 
quarks in the intermediate states, like the diagram in Fig. 3b. Formally, all such 
terms can be considered as the higher-order corrections to the kernel in the right- 
hand side of the equation (|3.1ip , 

G P (p\q) - G P (p\q) = GP( P \q)+\GP(p\q)+X 2 GP(p\q) + --- , (6.1) 

where Gp\p\q) = Gp(p\q) is the original kernel (|3.12j) . All terms here are expected 
to be singular at p = ±q, the singularities being of two types: the delta-function 
singularities coming from the disconnected diagrams (like the one shown in Fig. 3c), 
and the second-order poles at p ± q = reflecting the long-range character of the 
interaction in (13. 3p . Let us discuss first the role of these singular terms. 
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The delta-function terms were previously discussed in the Section 4. Their role 
is to modify the quark self-energy u(p) — > e(p), and to renormalize the quark mass. 
The first part is already taken into account in the left-hand side of the equation 
(|3,8p , and since the associated one-particle reducible diagrams disappear in the limit 
P = oo, this part plays no role in the final form of the Bethe-Salpeter equation (|4.8H . 
The second part - the quark-mass corrections due to the one-particle irreducible 
diagrams - can be taken into account if one replaces the parameter m in the Eq. (l4,8p 
by the full "dressed" quark mass m q . The latter is defined perturbatively, as the 
sum of all one-particle irreducible self-energy diagrams. By this definition the quark 
mass is given as the power series 

m\ = m 2 (1 + a 2 A 2 + a 3 A 3 + • • • ) . (6.2) 

In principle, the coefficients can be determined through well defined calculations 
within the perturbation theory. However, the problem is rather involved, and at the 
moment only the first correction term in (j6.2[) is known exactly [22], 

a 2 = 0.071010809... . (6.3) 



This and the higher order corrections in (|6.2p obviously modify the weak coupling 
expansions of the meson masses M n . Thus, in the low-energy expansions of the 
actual meson masses 

M 2 °° 

= 1 + ]T^ (6.4) 
4m z ^-^ 

k=l 

the coefficients ^ coincide with the corresponding coefficients e/% in (|5.13|) at k < 6, 
but starting from they receive additional contributions due to the quark mass 
corrections. In particular, at the order i 6 we have 

= £6 + 0-2 , (6.5) 

while the next term ~ t 7 receives no correction, so that we still have /jj = 0. Yet 
higher orders are affected by these as well as by other types of the multi-quark 
corrections, as we explain below. Our present knowledge about the "dressed" quark 
mass m q as the function of A is very limited. No higher order coefficients 03, 04, . . . 
in (|6.2p are known even approximately. Moreover, there are no reason to expect the 
perturbative expansion (|6.2p to converge. On a more fundamental level, no clear 
definition of the dressed quark mass beyond the perturbation theory is available. 
Therefore, apart from the weak coupling expansions, the quark mass m q can be 
treated presently as a phenomenological parameter at best. Taking this point of view 
we will observe in the Section 8 that the overall radiative correction 5m = m q — m 
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is probably very small (no more then ~ 2 — 3% of the lightest meson mass Mi) at 
all positive A. 

Next, the even order terms in (|6.ip are expected to have the second order poles 
at p± q = 0. The pole terms are interpreted as the corrections to the string tension: 
the "bare" string tension /o = 2a h in the r.h.s. of the BS equation (14. 8|) is effectively 
replaced by the "dressed" one / = f{h). To make it precise, we write 

foG P (p\q)=fGp(p\q) (6.6) 

where by definition the kernel Gp is normalized in such a way that its second order 
poles have the residues 1, as in f)3. 14|) : we can write 

G P (p\q) = G P (p\q) + AG p Cg \p\q) (6.7) 

where Gp(p\q) is the original kernel (13.121) which has the standard pole terms ()3.14j) . 
and the last term is regular at p = db q |zj. In this way the coefficient / appears as 
a power series in A 2 , 

/ = /o(l + c 2 A 2 + ---)- (6.8) 

It is possible to argue that the coefficients C2k here are related in a simple way to 
the coefficients <?2fc+i of the weak field expansion 

% - -L log m 2 = - \ A + ~ 92 A 2 + h A 3 + h A 4 + . . . , (6.9) 
m z 87T 2 

of the IFT vacuum energy density JT vac 0. Let us write T(h) for the analytic function 
of h which coincides with J- VSkC at positive real h. Then the analytic continuation 
of T(h) to the negative real values of h describes the energy density Tmeta of the 
"false" vacuum - the Lorentz invariant metastable state having the magnetization 
opposite to the external field h (e.g. (o"(x)) me ta < at h > 0). It is natural to 
interpret the difference 

^meta — ^vac (6.10) 

as the effective string tension. Indeed, when the two quarks constituting a meson 
are widely separated (the configurations responsible for the singular terms in (|6.7p ). 



16 Although regular at p ± q = 0, the term AGp BS \p\q) certainly has other singularities; even 
when p and q are real one expects to have branch cuts associated with the multi-particle thresholds. 

17 The vacuum energy density JT vac is of course the same as the (singular part of) specific free 
energy of the near-critical Ising model. 
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the spatial domain between the quarks is effectively filled with the "false" vacuum. 
It follows that the coefficients in (16.8P are expressed through g n as 



C2k = -2<72fc+i (6- 11 ) 



This relation provides several terms in (|6.8p . since many coefficients in (16.91) are 
known 0. In particular 

c 2 = -0.003889... (6.12) 

The correction terms in ()6.8p also make contributions to the low energy expansion 
(16. 4p . The first term affected is p% t 8 , 

, 2 C2 + ° 2 fa i q\ 

Ps = £8 H g £2 • (6.13) 

Corrections to the terms ~ i 9 and higher involve additional contributions from the 
regular term in (|6.7p . see below. 

It is tempting to take the difference (16.10P as the nonperturbative definition of 
the effective string tension /. It cannot be valid literally though. The function 
J-{h) is known to have a branch cut along the negative part of the real axis [23], 
hence the analytic continuation defining -F me ta in fact returns complex values. The 
imaginary part is interpreted in terms of the tunneling decay probability [24] and, 
being exponentially small in h, it is invisible in the perturbative expansion (|6,8p . In 
fact, both .Fvac and J r me ta. are known numerically with high precision, see Ref. [4]. 
In Fig. 6 we show a plot of the difference (|6.10p (measured in the units of |/i|i5) as 
the function of the scaling parameter (|1.2p . One can see that the imaginary part of 
(|6.10p is much smaller then the real part at all positive rj (moreover, it essentially 
vanishes at rj > 0.8), therefore taking the real part seems to provide a reasonable 
definition for the effective string tension /. Note that the real part differs very little 
from the bare string tension /o at all but very small (positive) rj. The difference is 
significant at ry < 0.8, and becomes drastic at rj —* 0. Indeed, at small 77 the bare 
tension tends to zero in a singular way, /o ~ rj s , while the difference (|6.10p has finite 
limit at rj = 0, and is analytic in some domain around this point. In fact, a number 
of terms of its Taylor expansion around r] = can be extracted from the results 
of [4]; for instance, for the real part / re of (16.10P we have 

/ re = \h\% [pv + Pirj + p 2 ri 2 + ■■■) , (6.14) 

where the coefficients are expressed through defined in Ref. [4] (see Table 3 there) 
as Pk = ( cos 8n ^ 5 ~ 2 ' ) — l) 3>fc , for instance 

po = 2.3692934..., p x = 0.3521342... , p 2 = , etc (6.15) 

18 The related coefficients G„ = (2s) n g n are collected in Table 2 of the Ref. [4] . 
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Figure 6: Plot of the difference (|6.10p measured in the units of \h\ 'A , vs the scaling param- 
eter r\. The solid line and the "fat" dashed line represent the real and the imaginary parts 
/ re = 5fte [^"mcta - ^"vac] and / im = 3m (^mcta - .^vac), respectively. The imaginary part is 
shifted upward by 1.6 to make the Figure more compact. The dotted line represents the 
bare string tension fo = (2sr}s) \h\~ , and the fine dashed line shows the effect of the first 
correction to fo explicitly written in Eq. (|6.8[) . 

According to the above analysis, the singular parts of the multi-quark corrections 
in (|6.ip modify the weak coupling expansions in a relatively trivial way, through the 
renormalizations of the parameters m and /o entering the BS equation (|4.8p . which 
are to be replaced by the "dressed" parameters m q and /. Thus, the singular terms 
are taken into account by substituting 

m — > m q , A — > f /mi (6.16) 

in the low-energy expansion (|5.16p . and the same applies to the semiclassical ex- 
pansion (I5.26|) as well. Remaining multi-quark corrections are associated with the 
regular term in (|6.7p . It is not difficult to show that this kind of corrections affects 
the low-energy expansion (16. 4p only at the order i 9 and higher. Indeed, the term 
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AG p 9 (p\q) is regular at sufficiently small p and q, and it is an odd function of each 
of this variables. Therefore its leading low-energy behavior is 

AG { p g \p\q) -» A{P) pq + --- . (6.17) 

Since this term is there entirely due to the higher-order corrections, the coefficient 
A(P) here is ~ A. Recall also that the low-energy domain corresponds to p, q ~ t. It 
follows that the correction to the BS equation due to this regular term is ~ t 9 = A 3 . 
Similarly, one can check that the regular term leads to the corrections to the semi- 
classical quantization condition (|5,26p at the order ~ A 3 and higher. Therefore 
developing the weak coupling expansions (both the low-energy and the semiclassi- 
cal ones) of the actual meson masses at the order A 3 and beyond requires explicit 
evaluation of the contributions from the "regular" term AGp Cg ^ in (16, 7h . The lat- 
ter constitutes a separate and rather involved problem which we do not attempt to 
handle here. Note that the correction term ~ A 3 to the quark mass in the Eq. (|6,2p . 
whose coefficient 03 is also presently unknown, contributes at the same order. To 
summarize, existing data are sufficient to determine the coefficients of the weak 
coupling expansions of M n all the way up to, but not including, the order A 3 . Eval- 
uating further terms requires additional ingredients, most importantly developing a 
perturbative expansion of the term AGp Gg ^ and the higher-order corrections to the 
quark mass. 

The above analysis, in particular the low-energy expansion, apply to stable 
mesons, with the masses M n well below the stability threshold 2M±, twice the mass 
M\ of the lightest meson. As given meson mass approaches to the threshold, one 
expects the threshold singularities of the last term in (|6.T[) to play more and more 
imoportant role. As the result, when close to the threshold the actual mass M n is 
expected to deviate significantly from its BS approximation M n . Furthermore, when 
the meson mass exceeds the threshold the actual mesons turn into resonances (see 
footnote 6 on page 4), and their masses M n develop imaginary parts representing 
the decay probabilities. This effect is certainly outside the scope of the two-quark 
approximation: the BS equation (|4.8p describes an infinite tower of stable mesons. 
And even including the multi-quark states perturbatively, as was discussed above, 
does not seem to provide a consistent way of treating the mesons above the stability 
threshold. Let us mention interesting attempt to treat the issue of the unstable 
mesons on the perturbative level in Ref. [13]. However, we believe that systematic 
treatment of the decay problem must involve at least some understanding of the 
meson scattering states, the problem which we hope to address in the future. 
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7 Numerical solution 



As we will argue in the next Section, the usefulness of the Bethe-Salpeter equation 
goes well beyond the derivation of the weak-coupling expansions (|5.16p . (|5.26p . 
Therefore it is of interest to study the eigenvalues M 2 of the BS equations (|4.1ip 
at all values of A, real and complex 0. The equation (|4.8p is unlikely to admit 
an analytic solution, but its numerical solution is not difficult to develop. The 
only feature which makes the numerical approach not exactly straightforward is the 
presence of the singularities of the kernel (14. 9p at 9' = ±6. This is one of many 
features that the equation (|4,8p shares with the mass spectrum equation in the 't 
Hooft's model [14]. Several methods have been developed for numerical solution of 
the latter (see [14,25]). We found it convenient to use yet another approach which 
is based on the Fourier-transformed version of the equation (I4.8p . We introduce the 
rapidity Fourier-transform 

/oo r-oo j 

il>(y)e- iv6 dv , ip{y)= I m(6)e wd —. (7.1) 
-oo J— oo 27T 

In terms of ip{v) the norm (|4.1U|) becomes 

1 f 00 

||^|| 2 = -/ il)*{v)K(v-v')ij)(v')dvdv' , (7.2) 



4 



oo 



where 



foo fj- v ® dQ 

K(v) = 1 ^ ^ = ^ffiW ' (7,3) 



and then the equation (14.81) transforms to 

8 ( m 2 + f o v tanh — ) i)(v) — — — — / T ib(v') dv 

V 2 / 2 coshf C o S h^ ^ V ; 

POD 

M 2 / [K(y - v) - K{v + v')} ij){v') dv 1 (7.4) 



Now the kernel is regular at any real arguments, and the equation (|7.4|) admits 
numerical solution through straightforward discretization of the variable v. Before 



19 There are good reasons to expect that the BS masses M n , taken as the functions of complex A 
(or rather of the complex variable 77 = 1/VX), imitate much of the analytic properties of the actual 
masses M n in the complex 77- plane. We will say more on this subject at the end of this Section, and 
again in Sect. 9. In our view, gaining any insight into the analytic properties of M n (rf) (structure 
of its Riemann surface, singularities, etc) is of central importance for overall understanding of IFT. 
This problem was one of primary motivations of our work. 
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turning to the numerics, let us say few words about the general properties of the BS 
equation written in this form. 

It is not difficult to show that generic solution ip{v) has poles at the values of u 
which solve the equation 

1 + A^tanh — = 0. (7.5) 
2 

where again A = /o /m 2 . At real positive A all solutions are purely imaginary, of the 
form ± i k, k being a positive root of the equation 

A k tan — = 1 . (7.6) 

Let kq be the lowest of such roots; at real positive A it lays in the interval [0,1]. 
The associated pole of ip(u) controls the large - 9 behavior of the rapidity-space wave 
function, 

V{6) -> r sign(#) e~ KoW as 1 | oo . (7.7) 

where r is a (normalization-dependent) constant. Note that at large A this pole 
approaches the real axis, since 

nr 

K o —> \ — r a s A — > oo . (7-8) 
V vrA 

At A = oo every eigenfunction ^ n [y) has a simple pole at v = 0; correspondingly, 
the associated rapidity wave function ^f n (6) tends to a constant r n at large \6\. One 
consequence of this phenomenon is the nature of the expansion around the point 
A = oo. This point appears to be a square-root branching point of the eigenvalues 
M 2 taken as the functions of A (see remark at the end of this Section). Therefore 
the large - A expansions of M 2 are of the form 

^| = a + r« Vx + + ... + y« ^ + ... , (7.9) 

(k) 

where Yn are constants, and the series converge in finite domains. It is straight- 
forward to derive the following expression for the coefficients at the terms ~ \/A, 

(7.10) 



y« 1 



2 V2vr 



Ir I 2 

'77, 



l*n| 



A=oo 



where r n are the constants in the asymptotics (17. 7p of the wave functions 
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n 


A = 0.01 


A = 0.10 


A = 1.00 


A = 10.0 


A = 100 


1 


1 


110840242 


1.550883388 


4.185654285 


24.29578962 


205.1229687 


2 


1 


196814699 


2.024246733 


7.637285443 


56.34010255 


523.5675928 


3 


1 


269079423 


2.446855032 


10.96393850 


88.01258645 


838.4720639 


4 


1 


334331298 


2.845183371 


14.23862629 


119.5937889 


1152.960249 


5 


1 


395100728 


3.228787972 


17.48445168 


151.1324905 


1467.318135 


6 


1 


452678928 


3.602360894 


20.71182429 


182.6463568 


1781.617148 


7 


1 


507833369 


3.968598364 


23.92635268 


214.1438592 


2095.883687 


8 


1 


561065333 


4.329208970 


27.13141319 


245.6299865 


2410.129928 


9 


1 


612724086 


4.685350869 


30.32919948 


277.1195629 


2724.362364 


10 


1 


663064726 


5.037848886 


33.52121875 


308.8383150 


3038.584826 



Table 1: Ten lowest eigenvalues of the BS equation (|4. at different values of A obtained 
by numerical solution of the Eq. (|7.4[) . The ratios M 2 /4m 2 are presented. 

As was already said, numerical solution of (|7.4jl is obtained by discretization of 
the variable v. In Table 1 we present few lowest eigenvalues for some values of 
A. At substantially smaller A (A < 0.001) the numerical results are indistinguishable 
(within the given accuracy) from the first eight terms of the low-energy expansion 
(|5.16p . On the other side, at large A, the expansions (|T.9[) apply. Numerical estimates 
of the three leading coefficients Y for few lowest mesons are given in the Table 2. It is 
interesting to observe that as n grows the leading coefficients Y^ quickly approach 
the following simple asymptotic form 

y n (0) -» vr (n - 3/8) as n -» oo . (7.11) 

In the same limit Yn tend to a constant value 1.2533..., while Y^ increases loga- 
rithmically as ^ log(n — 3/8) + Const, with Const ~ 1.209. 

Numerical diagonalization can be applied also at complex A. The procedure is 
straightforward as long as |argA| < tt. But at argA = ±7r the poles v = ±ikq 
surface at the real i/-axis, and if one continues to the next sheet |arg A | > tt these 
poles break through the real axis, leaving behind the residue terms. Thus, at the 
second sheet one has to consider the equation which differs from the Eq. (|7.4p by this 
residue term. Continuing around the point A = oo second time generates another 
residue term which cancels exactly the first one. It follows that A = oo is the 
square-root branching point of the solution, as was stated in the Eq. (l7,9j) above. 
Furthermore, one can argue that the BS masses M n , taken as the functions of 
complex A exhibit infinitely many singularities on the second sheet of the A - plane. 
Indeed, as was explained above, when continuing to the second sheet argA > tt, 
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n 


yr(0) 


1 n 


1 n 


n 


1 n 


yr(l) 
1 n 


V (2) 
- 1 n 


1 


1.904347602 


1.3810744 


0.873 


6 


17.67032120 


1.25432903 


2.072 


2 


5.094667187 


1.2671521 


1.450 


7 


20.81220388 


1.25403891 


2.154 


3 


8.242239347 


1.2582803 


1.691 


8 


23.95398627 


1.25385735 


2.224 


4 


11.38577501 


1.2558373 


1.853 


9 


27.09571021 


1.25373628 


2.286 


5 


14.52825003 


1.2548349 


1.974 


10 


30.23739767 


1.25365155 


2.341 



Table 2: Three leading coefficients of the expansion (|7.9p for the first ten eigenvalues M^. 
The Eq. (|7.10p was used in computing Yn . 



the pole at v = zko sinks into the lower half-plane, where at certain values of A it 
collides with other poles of ipfa), associated with the higher roots of the equation 
(j7.6|) . It is possible to show that these singularities accumulate towards the point 
A = 0, making this point an essential singularity of all BS masses. This picture was 
behind our statement in Section 5 about the asymptotic nature of the weak coupling 
expansions. But detailed understanding of the analytic properties of the BS masses 
M n as the functions of complex A remains to be achieved, and we intend to come 
back to this problem elsewhere. 



8 Comparison to the TFFSA data 

By construction, the BS equation was designed to describe stable mesons at suffi- 
ciently small A, where the idea of the mesons as predominantly two-quark constructs 
is well justified. However, it turns out that this approximation describes the mass 
spectrum of stable mesons remarkably well for all but very large values of A. More- 
over, when some multi-quark effects (notably the renormalizations of the string 
tension) are taken into account, the approximation works very well at all real pos- 
itive A. In this Section we compare the results derived from the BS equation with 
the data on the mass spectrum obtained directly from the IFT using the TFFSA [4]. 
It is convenient in this context to discuss in terms of the scaling parameter (|1.2p 
related to A as 



A 



2a h 



2s 

~. If ' 



.1) 



where s is the constant defined in (12. 2p . In the following discussion we still reserve 
the notation M n (or M n {r]) when the dependence on r\ is to be emphasized) for 
the actual masses of the particles in the IFT, and use the notation M n for the BS 



31 



masses, obtained from the equation (|4.8p . In fact, we will often write M n (m, /o) to 
signify their dependence on the parameters in the Eq, (|4.8j) , 

We do not describe here the TFFSA in any details, which are presented in 
Ref. [4] . The strategy for extraction the particle masses is not much different from 
that commonly used in the TCSA [12] . We have employed additional improvements 
(mostly in the way the finite-size effects are taken into account) which will be de- 
scribed in due details in a separate paper. We have used this technique to determine 
the actual masses M n of stable particles, at positive rj ranging from to 24. With the 
truncated space involving approximately 1800 lowest levels we were able to achieve 
typical accuracy of five significant digits for all those values of rj, except for when 
the mass M n is very close to the stability threshold 2 Mi, where the accuracy may 
be somewhat lower. The results were shown in Fig.l for rj £ [— 1 : 4] (at greater 
values of r\ deviations from the weak coupling approximations (|5.16p or (15. 26ft would 
not be visible). Beyond the threshold 2M\ the particles become resonance states, 
their masses M n acquiring imaginary parts[^|. At rj = most of the resonances 
disappear, while some become stable particles again, in accord with the spectrum 
of the integrable IFT at ij = 0. The masses at r] = are known exactly [5,26], and 
TFFSA reproduces them within the stated accuracy. For all we know, the integrable 
point T) = is analytic, in particular the masses M n = M n {rj) admit power series 
expansions 

M n (rj) = + M« n + t] 2 + ... . (8.2) 

The slopes for few lowest n are also known exactly, through the form-factor 

perturbation theory [6], and our TFFSA data is in full accord with the exact results 
(for previous numerical verification of the exact slopes using the lattice model see 
Ref. [27]). These numbers, as well as numerical estimates of the next term coefficients 
for the lowest mesons are collected in the Table 3. 

As expected, at large rj the masses from TFFSA agree very well with the low- 
energy expansions (|6.4p . This is demonstrated in Fig. 7, where the differences M n {rf)— 
M n (m, /o) (with r\ = (2sm 2 //o) 8 / 15 ) for the three lowest mesons are shown at large 
7] > 10. In this domain the BS masses M n (m,/o) are accurately described by the 
low-energy expansions (j5. 16H . According to the arguments in Sect. 6, the most part 
of the differences must be attributable to the quark mass renormalization (|6.2p . 

20 It is possible to extract the masses and the widths of the resonances from the TFFSA data, 
and we have preliminary results for few lowest resonances. However, so far precision of our data is 
far from being satisfactory, especially when the resonances are wide, and we refrain from reporting 
this data here. 

21 Shifting away from the integrable point renders unstable all but the three lowest particles; 
correspondingly, the coefficients (and higher) with n > 3 generally have nonzero imaginary 

parts [6, 16]. 
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n 








Mr 




m ] 


1 


4.404908579981 


1.295045 


0.2003 


4.248274145 


1.316490 


0.2815 


2 


7.127291799746 


1.115886 


0.2072 


6.948600321 


1.077762 


0.3726 


3 


8.761556059774 


1.953268 


-1.5110 


8.838163544 


1.095065 


0.3471 


4 


10.59322004129 


1.484334 




10.387731645 


1.144113 


0.3212 


5 


13.02221009790 






11.733999773 


1.253651 





Table 3: The coefficients of the Taylor expansions (|8.2j) for four lowest mesons. The leading 
coefficients are exact masses of the integrable theory at r\ = [5,26]. The exact slopes 
Mn are taken from Ref. [6], while the coefficients are estimated from the TFFSA 

data. The BS coefficients M^ 1 are computed using (|7.9p with A = //m 2 , and the expansion 
(l6d!) for /. 

since its effect on the meson masses is ~ A 2 , while other multi-quark corrections 
are of higher orders. At large r\ the deviations look consistent with the asymp- 
totic form (02 A 3 ) m = (4 02 s 2 l^l 15 ) r? -11 / 4 which corresponds to the leading quark 
mass correction (see Eq. (|6.5p ). The remaining mismatch visible in Fig. 7 is likely 
to come from the unaccounted (and presently unknown) higher-order corrections to 
the quark mass, and perhaps from the multi-quark corrections to the regular part 
of the effective kernel discussed at the end of the Sect. 6. 

It turns out that even if all renormalization effects are disregarded, the BS masses 
M n (m, /o) approximate the actual masses M n very well at all but very small r]. As 
is seen in the plot in Fig. 8, the BS masses M n (m, /o) significantly deviate from M n 
only at 77 < 0.5. Note that this is roughly the same region where the higher-order 
corrections to the string tension / become significant, see Fig. 6. The dramatic drop 
of the BS masses M n (m, /o) at rj — > visible in the plot is certainly related to 
the singular behavior of /o in this limit. Therefore it is reasonable to expect that 
replacing /q with renormalized string tension / would take care of the most part 
of the mismatch at small rj. As was already said, it is not completely clear how to 
determine the string tension / beyond the perturbation theory, but it is unlikely we 
can be off too far if we take, say, the real part of the difference (|6.10p . i.e. set 

/ = /re = ^6 (^meta — -^vac) • (8.3) 

In particular, we recall that / re has finite limit at rj = 0, and in fact analytic around 
this point. As the result, the BS masses M n (m, f rc ) are also analytic at 77 = 0, 
and their expansions around this point follow from (|7.9p (where this time one sets 
A = / re /m 2 ), and <KWf . We show some coefficients of these expansions in Table 3, 
to compare with the coefficients in (18. 2p . To get further feeling of the matter, we 
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Figure 7: Plots of the differences M n — M„(m,/o) (in the units of vs the scaling 

parameter r\. The solid line is the curve (4 02 s 2 ) 77 -11 / 4 representing the correction 77102 r 6 
from the term 02 in Eq. (|6.5p . 

have computed numerically few lowest BS masses M n using / re instead of /o in (|4.8|) ; 
the results (also shown in Fig. 8) indeed come rather close to the actual masses at 
all positive 77. It is tempting to attribute the remaining mismatch of the lowest BS 
masses M n (m, / re ) (clearly visible in Fig. 8 for n = 1,2, 3) to the effect of the quark 
mass renormalization, so far disregarded. Indeed, if the correction to the quark mass 
is positive (as the A 2 term in (|6.2p is), its effect would be to shift the meson masses 
upward. Unfortunately, nothing is known about the higher-order corrections to m q 
(other then the term (|6.3p . that is), let alone anything beyond the perturbative 
expansion. Just adding the known A 2 term, while noticeably improving the BS 
masses at sufficiently large r\ (as was seen in Fig. 7), leads to substantial overestimate 
of the actual masses at small rj (Fig. 8). This suggests that the properly defined m q 
remains close to the bare quark mass m at all positive rj, being perhaps no greater 
then (0.1) |/i| 8 / 15 at rj = 0. 

Since the two-quark approximation can not account for the possible instability 
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Figure 8: Plots of few lowest meson masses (in the units of | | "ti" ) , vs the parameter r\. The 
bullets • represent the actual masses (exact values at rj = 0, and the TFFSA data elsewhere. 
The uncertainties in the latter data are smaller then the size of the bullets). The fine dotted 
line is the stability threshold 2M\. The dashed lines are the masses M n (m, /o) from the 
BS equation with no renormalization effects taken into account. The solid lines are plots of 
M„(m, / re )- Effect of the first correction term in (|6.2| is shown by the dotted lines. The 
circles o indicate position of the higher thresholds M\ + Mi and Mi + M3 at r\ = 0. 

of the heavier mesons with respect to decays into the lighter ones, the BS masses 
M n remain real at all positive 77. This is unlike the actual mesons which turn 
to resonances as soon as their masses cross the stability threshold 2M±. In fact the 
behavior of the actual masses M n {rj) near the threshold is rather subtle (see footnote 6 
on page 4 for comment); this complex threshold behavior of the masses M^{r]) and 
M^iyj) is even visible in Fig. 8. Of course this subtleties are not captured by the BS 
equation (14, 8j) - the BS masses M n go right through, without any anomalies at the 
threshold. This deficiency is the reason why in the Table 3, while is rather close 
to Mg ^ , the slope and the curvature are so much off - the point r\ = 

is too close to the stability threshold 7/3 = —0.136(1) of the third particle. Let us 
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recall that it is when the meson mass approaches the threshold that the multiquark 
corrections to the "regular" part of the BS kernel, Eq, (|6.7j) . are expected to play a 
prominent role, and it is plausible that including such multiquark corrections would 
help to understand better the threshold behavior. 

In this paper we do not discuss the resonance part of the spectrum, putting this 
task away for another time. However let us make few remarks here. Although the BS 
equation (|4.8f) can not account for the instability of the mesons above the threshold 
(and thus disregards the imaginary parts of their masses), it likely captures the most 
of the real parts of the resonance masses. This is not inconsistent with the pattern 
readily observed in Fig. 8. For example, the BS masses M^m, f re ) and M§(m, f re ) 
at rj —* come reasonably close to the exact masses of the fourth and the fifth of the 
stable particles of the rj = theory, suggesting that analytic continuations of the 
actual masses M±(rj) and M§(rj) to the point rj = reproduce those exact masses. 
This rises a question about the fate of the mass M^{rj) at rj — > 0. Now, observe that 
the BS mass M^{m, / re ) at rj — > lands far from (this is cleaarly seen in the 

last line of Table 3) but rather close to m[ + , the sum of two lowest masses 
of the rj = theory. The most economic scenario then implies that for the actual 
masses we have Ms(0) — M\(0) — M2(0) = 0, and at rj « 1 the state with mass 
M§(rj) can be interpreted as a weakly coupled bound state of the two particles M\ 
and M.2- Of course it would be interesting to elaborate quantitative details of such 
scenario (and its versions for the higher resonances); again, we hope to address this 
problem in the future. 

9 Discussion 

In this paper we refine technique which allows for quantitative understanding of the 
mass spectrum of the IFT, in the low-T domain, in terms of "mesons" consisting 
predominantly of two quarks bound by a confining interaction. Thus we develop 
further the idea originally due to Wu and McCoy [3]. Systematic approach based 
on the idea of the mesons being essentially the two-quark constructs leads to the 
Bethe-Salpeter equation, Eq. (j4.8j) . yielding a discrete spectrum of the meson masses. 
The Bethe-Salpeter equation is itself an approximation since it essentially ignores 
multi-quark components of the bound states. The approximation becomes exact in 
the weak field limit rj —* oo, but we show that it reproduces the masses of stable 
particles with reasonable accuracy at all positive rj. 

There is much room for further development. First, it seems to be possible 
(albeit involved) to take into account the multi-quark corrections perturbatively. 
This problem was briefly discussed in Section 6. It would be especially important 
to understand the role of the multiquark corrections to the "regular" part of the BS 
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kernel, the term AGp in Eq. (|6.7p . That would help to get at least some insight 
into the way the meson masses behave near the stability threshold. Understanding 
of the above-threshold (i.e. resonance) part of the spectrum is another open problem 
(see however [13]). 

As was mentioned in the Introduction, the main motivation of this work origi- 
nated from our attempts to understand analytic properties of the IFT masses M n (rj) 
as the functions of complex rj. There are reasons to believe that at least some singu- 
larities of these functions are mimicked by the BS masses M n (m,f), considered as 
the functions of complex variable fj = mj^f] ', and understanding analytic properties 
of the BS masses is a much more manageable problem. For instance, one type of 
singularities - those appearing due to collisions of the roots of the equation (|7.5|) - 
was mentioned at the end of Sect. 7. These are square-root branching points, located 
somewhere at |arg7/| > 37r/4 (We plan to discuss these singularities in more details 
elsewhere). It is likely that the actual IFT masses M n {rf) have similar singularities 
at | arg 77 1 > ll7r/15. Another question of much interest concerns locations of possi- 
ble zeros of the BS masses M n on the Riemann surface of fj. At least some of such 
zeros may signal divergence of the correlation length in the full-fledged IFT, where 
one would expect to observe associated critical points. In this manner, vanishing 
of the lowest meson mass Mi (77) at r/ = (2.4295..) e ±l ^5 n is related to the Yang-Lee 
edge singularity [4]. It is certainly interesting to contemplate the possibility that 
the Yang-Lee edge singularity is but one of many critical points one can bump into 
when wandering on the Riemann surface of r/ in the IFT. To some extent, the present 
paper can be regarded as a preparation to this analysis, and we hope to address this 
family of problems in future work. 

Finally, some of the techniques developed here, in particular the weak coupling 
expansions in Sect. 5, can be useful in other 2D theories which exhibit confinement 
of quarks, notably the 2D model of multicolor QCD (the 't Hooft model). More 
generally, there is a large class of systems having degeneracy of the ground state, 
which are integrable (with factorizable S-matrix of kinks), but not free. Typical 
perturbations of such theories lift the ground state degeneracy and creates confining 
force for the kinks. It seems important to extend the approach developed here to 
such systems. 
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Appendix 



A Infinite-momentum limit of Eg. (13. Ill ) 

Here we derive the P — * oo limit of the equation (|3.1ip . We assume that e(p) has 
exact asymptotic behavior f)4.6f) . and that AP(P) has the form (|4.5j) . in particular 



M 2 

AE = \P\ + — + 0(\P\- 3 ) . 



(A.l) 



Both sides of (|3.1ip can be expanded in inverse powers of P. Assume that \Pp(p) is 
normalized in such a way that at \p\ < P/2 it has finite nonzero limit as P —* +00. 
Under such normalization, ^p(p) with |p| > P/2 tends to zero as P —* +00. More 
precisely, if it = 2p/P is fixed, 



fO(l) 

\0(P-2) 



at 
at 



|n | < 1 
lul > 1 



as P 



oo . 



(A.2) 



Indeed, consider the left-hand side of the equation (|3.1ip . The coefficient in front of 
^>p(p) has different P — * oo behavior depending on whether \u\ is smaller or greater 
then 1, namely 



l.h.s. of (3.11) 
where 



p-1 r 2m? 
ll-u 

P [ 



— — —y | $(«) at I u I < 1 
u I - 1 I $(V) at I u \ > 1 



= lim tfp(uP/2) . 

P^oo 



as P -> oo , (A.3) 



(A.4) 



On the other hand, the kernel Gp(p\q) decays as 1/P 2 at all p, q ~ P, and hence 
the right-hand side of (|3.1ip is ~ P -1 at all u. As the consequence, the limiting 
wave-function $>(u), Eq. ()A.4p . must vanish outside the interval \u\ < 1, 



$(it) = at \u | > 1 . 
Comparing the 1/P terms in both sides of (|3.1ip . we obtain the equation 



m 



1 



M 2 

4 



where 



F{u\v) 



^/(l-u 2 )(l-v 2 ) 



1 — 

(u — -u) 2 



1 + uv uv 
(u + v) 2 4 



(A.5) 



(A.6) 



(A.7) 
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Changing variables u = tanh#, v = tanh#' one arrives at (|4.8p with *S>(9) = 
3>(tanh<9). 

Note that the equation (|A.6P is very similar to the Bethe-Salpeter equation 
determining the mass spectrum of mesons in the 't Hooft's model [14]. When the 
two quarks have equal masses the latter differs from (|A.6|) only in the detailed form 
of the kernel, which in that case is somewhat simpler, 

-F'tHooftHv) = 7 — — (A.8) 
[u — V) z 

(also, the parameter /o is to be interpreted as the square of the gauge coupling 
constant). 



B Derivation of the Eq. (15.51) 

Here we evaluate the action of the operator H — M 2 (with H defined in (14.120 ) on 
the wavefunction (|5.2p . It is useful to split the letter into two pieces, 



*o(0) = - , (B.l) 

with 

and evaluate the corresponding two terms separately. In fact we need only to evalu- 
ate \H — M 2 ] ^+(9); the desired result is then obtained by taking the part antisym- 
metric with respect to the reflection 9 — * —9. It will be also convenient to represent 
the operator H — M 2 as the sum of two terms 

[H - M 2 }^{9) = Am 2 cosh 2 9 [ Cw(d) + G^{9) ] , (B.3) 

where 

mo) = mm, with m = -^s(o) = i-^— ^, ( B . 4 ) 

and 

In the last equation we have introduced the notations 

V ; sinh 2 (0-0') J -oo cosh 2 9 KJ 2tt v ; 
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Let us first transform the integral 



poo 



xj a 

J — oc 

d_ 

d& 



01 



riff' 



(B.7) 



A 



dff 



sinh(6>-6") 2tt 



e *S(/J) 



d/3 



sinh(6" - p - iO) 



The order of integrations can be interchanged, and at fixed real /3 the contour of 
integration over 8' can be shifted downward, 8' = a — in/2; this leaves behind half 
of the residue at the pole 8' = 8, 



(Bl 



8 A ei s ^d(3 t d_ 
88 smh(8 + - iO) + dd 



00 dp 



A ei s ^da 



cosh(# — a) cosh(a + 0) 



In the first term we transfer the derivation over 8 to the integration variable (3, and 
then integrate by parts. The integral over a in the second term is evaluated in closed 
form. One finds 



(Bl 



1 



£1(0) ei s ^ dp A 
2 7_ 00 sinh(6> + P - iO) + tt 



d 8 + p 



dp sinh(# + P) 



dp , (Bi 



where the factor — fl(P) in the first term came from the derivative of S(P), see 
Eq. (|B.4|) . Note that when plugged into (|B.3|) . the first term combines nicely with 
the term £l(8) 1 $> + (8), to produce 



- n(p) eim 



i 



2 sinh(6» + p - iO) 



d/3. 



(B.9) 



Since the numerator in the integrand vanishes at /3 = —8, this point is no longer 
singular, and the shift — iO in the denominator is irrelevant. In fact, the integrand 
in (|B.9P can be reduced to 



M 2 



sinh(fl - p) 



8m 2 cosh 2 8 cosh 2 /? 



(B.10) 



Next, when applied to the second term in right-hand side of (|B.5|) involves 
the factor 



sinhfl dd 



cosh 2 8 2tt 



S(J3) 



, sinh((9 + /3 - iO) 
sinh/3 1 d p 



dp 



4 cosh 2 /? 2vr 9/3 cosh /3 



dp, (B.ll) 
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where the expression in the second line results from explicit evaluation of the integral 
over 9. 

Putting the Eq's. (|B~4|) . ([R5jl . (lR8l) . and ([Bill together, we find 



[ft + G]^+(#) 



M 2 sinh((9 - (3) 



+ 



8m 2 cosh 2 6> cosh 2 /3 
\ d f 9 + (3 1 sinhfl ivr/2 + /3 

0/3 Vsinh (# + /?) ~ 8 cosh 2 9 cosh/? 



7T 



ei 5 ^d/3. (B.12) 
Finally, isolating the odd part of this expression, one arrives at (I5.5p . 
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